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Abstract 

The broken symmetric phase of scalar models exhibits an infrared fixed point 
which is induced by the degenerate effective potential. The definition of the 
correlation length in the infrared regime enables us to determine the type of 
the phase transition in the model. It is shown that the massive sine-Gordon 
model exhibits a continuous, while the layered sine-Gordon model has an 
infinite order Kosterlitz-Thouless type phase transition. 



1. Introduction 

The renormalization group (RG) method provides a powerful tool to map 
the phase structure of field theoretical models [H, 0] • The phase space usually 
contains singularity points, which correspond to stationary solutions of the 
RG flow equations. The linearization of the RG equation in their vicinity 
results in the power-law scaling of certain physical quantities which defines 
the corresponding critical exponents j^]. In order to get the critical exponent 
v of the correlation length £ one usually should solve the stationary RG 
equations around some crossover fixed points Q. In the d-dimensional O(N) 
model the Wilson- Fisher (WF) fixed point while in 2d sine-Gordon (SG) 
model @, 0, 8 the Kosterlitz-Thouless (KT) fixed point @, El} plays the 



role of the crossover fixed point [111 ]. However there exists several models 
where the phase space either does not contain any crossover, saddle point 
like fixed point. These models also can possess non-trivial phase structure 
even including a broken symmetric phase. In these cases it is extremely 
difficult to determine the order of the appearing phase transitions. 
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Among many examples the massive SG (MSG) model has not any fixed 
points due to the trivial scaling of the mass. The model is a simp le gener- 
alization of the SG model with an auxiliary mass term 12J, [l3|, |l4], |l5j. The 



UV scaling of the MSG model predicts two phases with an infinite order 
phase transition, since one cannot distinguish the model from the original 
SG model there. Otherwise it is well known that the model as the bosonized 
version of the 2d quantum electrodynamics (QED 2 ) [l^] possesses a finite 
order phase transition. We note that the QED2 is one the most widely used 
toy model of self-interacting gauge theories since the model exhibits the soft 



mechanism of the confinement [17|, [12J , and it can be related to the chiral con- 
densate 18j. The phases of the bosonized version of QED 2 was first mapped 
with RG method in [13j, which successfully recovered the phase structure of 
the original fermionic model, moreover it was shown that the critical ratio 
u/M 2 in the IR regime exhibits a quantitative agreement with the value ob- 



tained by density matrix RG method [19|. The functional RG calculation 
was performed in the lowest order of the gradient expansion i.e. in the lo- 
cal potential approximation (LPA). This result was based on considering the 
upper harmonics of the periodic term in the MSG model, furthermore it was 
assumed that the RG flows run into singularity 13|, |20) . However the phase 



structure can be uncovered without taking into account the higher harmon- 



ics 



151 ]. although the inclusion of wavefunction renormalization disables us 



to map the scalar model onto QED 2 , due to the scaling of the parameter 



j3 [15J. The mass term of the MSG model on one hand breaks the periodic 
symmetry, and, on the other hand sets the RG evolution of all the couplings 
relevant below the mass scale. The broken periodic symmetry brings the KT 
type transition of the SG model into a second order one in the MSG model. 
It is believed [l9| that the model belongs to the 2d Ising universality class 
due to its Z(2) symmetry. Here we show that the exponents does not give 
the 2d Ising ones suggesting that the model may not belong to the 2d Ising 
universality class. 

The layered SG (LSG) model is an other generalization of the SG model, 
which contains N scalar fields interacting through a mass matrix 2l[. The 
model preserves the periodic symmetry of the original model. The LSG 
model is the bosonized version of the multi-flavor QED2, where the number 
of layers of the LSG model equals the number of flavors in QED 2 . The LSG 
model can also be used to describe the vortex behaviour of magnetically cou- 
pled layered superconductors, e.g. artificially produced superlattices, where 
the vortices interact through the transverse magnetic field [22J. In this treat- 



2 



ment one should assume that the Josephson coupling due to Cooper pair 
tunneling between the superconducting layers is suppressed in the model by 
the relatively large insulating layers. 

The LSG model also exhibits two phases as was demonstrated by using 
the functional RG method in an extended UV regime [21.]. It was also shown 
that the phases are separated by a critical parameter (3 C which depends sim- 
ilarly on the layer number iV as the critical temperature of fractional-flux 
in magnetically coupled high-temperature superconductors. The critical pa- 
rameter separates the symmetric and the broken symmetric phases of the 
model. It is assumed that there is a KT type phase transition between the 



phases [23], but an LPA investigation cannot account for that. It was shown 



in the case of the SG model 24 that one should include the next order of 



the gradient expansion, i.e. the wavefunction renormalization |25[ in order to 



unfold the KT-type crossover scaling regime [26|, |27|, [24J , therefore it seems 
unavoidable to go beyond LPA to establish the KT like character of the 
phase transition in the LSG model. However, similarly to the MSG model, 
RG flow equations for the LSG model have no crossover fixed point because 
of the trivial evolution of the interlayer coupling. We note that the due to 
the inclusion of the wavefunction renormalization the parameter (3 evolves, 
which makes impossible to map the LSG model onto the multi-flavor QED2. 

The RG flow equations of the MSG and the LSG models have no crossover 
fixed points disabling us to determine v in a usual way jij. However it 
was shown in several scalar model examples 11, 0] that the IR scaling of 
£ coincides with the one found around the crossover fixed point due to the 
global nature of the appearing; condensate in the broken symmetric phase 
28]. Our method in 11, 0] was capable of describing both the Ising-type 



and the KT-type phase transitions and scalings giving us a possibility to find 
the exponent v in the IR limit. There we define the correlation length as the 
reciprocal of that momentum scale k where the evolution stops in the broken 
symmetric phase due to the condensate induced degeneracy. We determine 
the type of the phase transition and the corresponding universality class 
according to the scaling of £ and the value of the exponent v. We emphasize 
that this method is the only possibility to determine the type of the phase 
transitions for the MSG and the LSG models directly. Furthermore the usual 
technique based on characterizing the phase transitions around a crossover 
scaling regime is not reliable, since the RG flow should collect all the quantum 
fluctuations during the evolution therefore one should follow the RG flows to 
the deep IR region. This fact makes the finding of any infinite order phase 
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transition rather accidental for example in the SG model around the KT type 
fixed point. Fortunately the crossover and the IR scaling regimes provide the 
same exponent v due to the global nature of the condensate which is closely 
related to the correlation length in the broken symmetric phase. 

For the MSG model we find a second order phase transition character- 
ized by the exponent v = 0.5. This suggests the possibility of exhibiting a 2d 
Ising-type phase transition there, which should give v — 1. We investigate 
the phase structure of the LSG model in a great detail. We show that the 



LPA approximation seems to show two phases in accordance with [17|, Il2 



However the inclusion of the wavefunction renormalization changes drasti- 
cally the IR behaviour of the model, and the broken symmetric phase melts 
into a single phase. We also determine the scaling of the correlation length 
which proves that the model possesses a KT type phase transition. 

The organization of this paper is the following. First the scaling in a finite 
order phase transition is investigated in the framework of the MSG model in 
Sect. [2j In Sect. [3] the LSG model is discussed. Finally, Sect. @] contains the 
summary of the results. 



2. Finite order phase transition 



The successive elimination of the quantum fluctuations is performed by 
means of the the Wetterich RG equation for the effective action [l[ 



* l m R 
T = -Tr 

2 R + T" 



(1) 



where ' = d/dip, ' = kdk and the trace Tr denotes the integration over all 
momenta. Eq. ([T]) has been solved over the functional subspace defined by 
the ansatz 



(2) 



where the potential V has the form 

V = -m 2 ip 2 + u n cos(rup), 



N v 



(3) 



n=l 



for the MSG model, with mass m = m^, the coupling u n = u n ^ and the wave- 
function renormalization Z = Z^i^p). We introduce the higher harmonics of 
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the periodic term of the potential, since these terms are generated during the 
RG evolution. The ansatz for the wavefunction renormalization Z contains 
a constant and a periodic term, i.e. 



Z = z + 2_. z n cos(nip). 



(4) 



n=l 



We note that, as opposed to the polynomial models, the field independent 
coupling zq has a non-trivial evolution without considering the upper har- 
monics z n , n = l...N z . The polynomial suppression has the form 



(5) 



with b > 1. The choice 6=1 coincides with the usual Callan-Symanzik 



scheme 29]. From the functional RG equation in Eq. (TjQ) one can deduce 

and z n - For shorthand we 



u, 



evolution equations for the couplings m 2 
introduce u = u\ and z = zq. The evolution of mass decouples from the 

where A is the UV 



ni 



other RG flow equations, and one obtains m 
cutoff. The further couplings scale according to the equations 



A- 



U, 



V n / RG 



-Z"G 2 + ( - d z% 2 + AZ'F'" ) G 3 



r'" 2 ( d p2 p + -/9 2 p\ + - d Z'p 2 T'"d p 2P 



+-,p 2 T" l2 d p 2P 2 G 5 
a 



G l 



(6) 



with G = l/(Zp 2 + R + V") and P = Zp 2 + R. We introduced the projection 
as V n = f cos(rup) / V for n > 1 and Vq = The momentum and the 

projection integrals have to be performed numerically. The scale k covers the 
momentum interval from the UV cutoff A to zero. Typically we set A = 1. 

The RG evolution becomes singular, i.e. the r.h.s. of the RG flow equation 
is infinite in Eq. (|6]) when 



Zp 2 + R + V" = 0. 



(7) 
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The non-trivial expectation value emerges from the average (ip) = contin- 
uously, therefore one can set <p = 0, giving 



k 2 -^2n 2 u n = 0, (8) 



ra=l 



where k = min P has been introduced. By using the polynomial suppression 
in Eq. ([5]) one arrives at 



P = b P ( ^n=0 Z n ^ (g) 



for 6 = 1 we have k = k. We introduce furthermore the normalized funda- 
mental mode u n = u n /k 2 . 

2.1. Infrared scaling 

We numerically calculated the flows of the coupling u in the presence of 
power-like suppression scheme. The mass scales trivially, i.e. fh ~ k" 1 . We 
considered N u = 5 and N z = 1 for the number of couplings. 

The coupling scales as u ~ k~ 2 below the mass scale in the IR region, 
independently on any initial UV microscopic parameters. However the model 



has two phases [13|, [15|, ll9(, separated by the critical ratio uo c /ml, where 
uq c is the critical IR limiting value of the coupling u, which separates the 
appearing two phases in the model. Since the mass m does not evolve, then 
m = rriQ. The phases can be distinguished by considerin g th e sensitivity of 
the effective potential to the microscopic parameter ua 



symmetric phase the effective potential depends on the UV value u a, while 
in the broken symmetric phase the effective potential is is insensitive to any 
parameters, i.e. it is (super) universal. 

Setting the value of the mass m fixed one has a broken symmetric phase 
when ma > «a c- Fig. [1] shows that one can also distinguish the two phases 
of the model according to the scaling of z. In the symmetric phase z is 
constant in the UV, scales as z ~ k~ Vco in the crossover and is constant 
too in the IR limit. The crossover scaling appears below the mass scale far 
from the deep IR regime. The anomalous dimension rj being characteristic 
for the divergence of the correlation function of the field variables, defined 
as rj = —dlogz/dlogk is not unique during the evolution, and its value is 
sensitive to the parameter choice of b, too. The final constant value of z in 
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the IR limit diverges as the initial parameters approach the critical ratio. 
In the broken symmetric phase we can find all the three regimes for z as 
in the symmetric case, but the IR regime shows its relevant scaling with a 
high value of anomalous dimension tjir. We also plotted the evolution of 
z as the function k. In the symmetric phase the flow is qualitatively same 
that was obtained for k. However the crossover scaling gives universal, i.e. b 
independent exponent fj co = 1 and in the broken symmetric phase the flows 
blow up in the IR, giving fjm = oo. This defines a dynamical momentum 
scale k c , where the RG evolution stops and the effective potential becomes 
degenerate. We got similar behaviour for other scalar models, too [HI, 0|. 
The reciprocal of k c is identified as the correlation length, i.e. £ = l/k c . The 
initial UV value of z Q \ can be identified with the temperature, therefore its 
distance from its critical value (for a fixed value of u\) provides the reduced 
temperature t. This makes us the possibility to determine the t dependence 
of the correlation length. In Fig. [2] one can see how £ depends on the critical 
temperature t. We performed the calculation for N u = 5 and N z = 1 and 
obtained that the correlation length scales as a power-law function giving 
straight lines in log-log scale according to 

f~r" (10) 

with v « 1/2. We repeated our calculations for the case N u = 1 and N z = 
and and for several values of b, and we got the same value for the exponent. 
This means that neither the upper harmonics nor the field dependent wave- 
function renormalization Z affects the value of v. This result is in accordance 
with the ones obtained for the SG model where we also got such sensitiveness 
of v to the upper harmonics. We also made calculation in the case of LPA, 
which also gave v = 1/2. 

Our results seems to contradict with results obtained for the 2d Ising 



exponent, where v — 1. In |19[ it was shown that QED 2 is in the 2d Ising 
universality class, and it is argued that the MSG model has a Z{2) symme- 
try which implies that the model is in the 2d Ising universality class, too. 
Although the MSG model has only a Z(2) symmetry, the RG flow equations 
for the couplings u n seems to preserve the periodicity. This may affect the 
simple Z(2) symmetry and may eventually give that the MSG model is not 
in the 2d Ising the universality class. 

These results demonstrate that one can get a clear picture from the type 
of the phase transition in the MSG model according to its IR scaling. Since 



7 



the model has no crossover fixed point this method is the only possibility to 
uncover the order of the phase transition there. 

Considering the RG evolution equations in Eq. © one can easily see that 
we have no "static" solution, which means that in mathematical sense we 
have no nontrivial fixed points in the MSG model. So far we investigated 
such models where the IR value of v was inherited from its crossover value 
11, if. In the MSG model it is not the case. This fact clearly shows the power 



of finding the scaling relations in the IR. We conclude that the IR scaling 
gives us an exclusive and correct picture for the scaling of the correlation 
length. 



3. Infinite order phase transition 

The LSG model can be described by the the corresponding effective action 



30 



r fc 



N 



71=1 



with the scalar fields of the layers (p n , n — 1 . . ., the wavefunction renormal- 
ization z, and the potential V of the form 



^ / N \ 2 N 

V = 2 J E<m + M E 

\rt=l / n=l 



COSLp n , 



(12) 



with J the interlayer coupling and a single coupling u for every layer. The 
elementary excitations with higher vorticity, which are coupled to the op- 
erators like cos(nyj m ), with n > 1 and m = 1...N are negligible in the 
determination of the thermodynamic properties of the model. Otherwise, 
the fundamental mode u is capable of mapping the phase structure of the 
model 2l|, |l4], |3l| . We consider only a field independent wavefunction renor- 



malization constant z, the field dependent ones do not affect the results in 
the case of periodic models. 

In order to find the flows we use the same Wetterich equation in Eq. (TT]) 
with the matrix 

T" = zp 2 \+V" (13) 



and V" 



5 2 V k /6<pi6<pj, i,j 



N and now the trace Tr denotes the 



summation over the internal index (over the components of the scalar field) 
and the integration over all momenta. Now we set b = 1 in the suppression 
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if the wavefunction renormalization z also evolves, since although it gives 
moderate convergence and needs high computational accuracy, it makes the 
evolution equations much simpler. The RG evolution equation in Eq. (JTJ 
gives differential equations for the dimensionful couplings z, u and J. The 
latter evolves trivially i.e. 

j = 0. (14) 

For the potential we have 

V^^G ] = ^l zpn + R m+v „ (15) 

where tr is a trace over the components of the scalar field and 11 denotes the 
iVth order unit matrix. The evolution of u can be picked up by the acting 
on Eq. (1151) with the projection operator 

V N = [ cos y2 2 ... cosy? at, (16) 

(2vr) JV J V1 ... VN 

to get the corresponding Fourier mode. The evolution of z becomes 



TV 



z = £ I / khl [G P V" [n) G P V" (n) (2pW-zG 2 p )]. 

n=l J P 

(17) 

for the choice b = 1, where V"^ = 5V" /5tp n . 

3.1. Local potential approximation 

We solved the RG evolution equation in Eqs. (TH)) . (fTo"|) and ffTT]) for 
the LSG model numerically. We determined the RG flows by taking into 
account only a single mode in the potential, but without taking any further 
approximations. The calculation is carried out for bi-layered systems, i.e. 
N = 2. This specification on one hand is necessary because the momentum 
integral and the integrals in the projection Vn has to be performed with no 
truncation in u in order to get the deep IR evolution. On the other hand the 
qualitative results are not expected to change by increasing the number of 
layers. 

We are interested in the LPA results because the bosonization technique 
makes the LSG model equivalent to the multi-flavor QED 2 . It is assumed 
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that the latter model has two phases. The mapping between the models 
can be performed by the parameter choice z = 1/4-7T. This requires the 
freezing of the wavefunction renormalization and confines our treatment to 



LPA approximation [31 



Acting on the projection V2 to the evolution equation for the potential 
in Eq. f lT5|) . one obtains 

^ _ f R(2P - m(cos0i + cos0 2 )) cos ^2 
Jp^z (P -u COS (fii){P -u cos ip 2 ) - J 2 

with P = zp 2 + R + J. The interlayer coupling J does not evolve, J = Ja- 
In the LPA treatment the evolution of z is also frozen. The extended UV 



scaling gives [2l|, [31 



1/8ttz / ^2 _|_ 2 J \ V 167r2 



" = "HaJ U+wJ ' (19) 

The coupling u = u/k 2 scales irrelevantly (relevantly) in the symmetric (bro- 
ken symmetric) phase if 1/z > 16tt (1/z < IQtt), respectively. We note, that 
generally, the critical frequency which separates the two phases of the model 
is [21 

z c = ^— ^. (20) 

In the broken symmetric phase the RG flow of the couplings can run into 
singularity if the numerical calculation is not precise enough. For example in 
the case of the SG model the exact treatment in the couplin g u is necessary 



to avoid the singularity for certain suppression schemes [24], |21| . Therefore 
we should include all the contributions of the coupling u exactly for the LSG 
model, too. 

Similarly to the MSG model we introduce k = min P, the dimensionless 
couplings u = u/k 2 and u = u/k 2 . The latter tends to 1 if the flow approaches 
the degeneracy. 

3.1.1. The phase structure 

The extended UV results suggest that the model has two phases. In 
Fig.[3] we plotted the trajectories of the phase space at a certain value of 
interlayer coupling J = 10~ 4 . We should admit that the phase space is 3- 
dimensional, but J scales trivially, therefore we always project the phase 
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space to the other axes. There is a symmetric phase at z < 1/167T. The 
flow of the coupling u is irrelevant, i.e. it tends to zero in the deep IR 
regime. Although the flows should span the momentum interval k G [A, 0] 
the evolution cannot reach its lower limit, because we should keep Ak/k small 
throughout the evolution. We note however that only the deep IR scaling 
of the couplings which provides us all the information about the effective 
potential of the model. The numerical results show that in the deep IR the 
coupling is sensitive to its initial value u\. There is a line of fixed points 
in the symmetric phase at z < 1/167T and u — 0, which is denoted by thick 
straight line in Fig.|3j 

When z > 1/167T, then the RG flow of the coupling u belongs to the the 
broken symmetric phase. Fig. [3] shows that if 1/87T > z > 1/167T then u tends 
to an IR fixed point denoted by a thick solid line there. We note that the 
concrete value of these fixed point is scheme dependent, as is demonstrated 
in Fig.Hl This behaviour is quite similar to that was obtained in the case 



of the SG model [11] . The singularity is avoidable by an appropriate choice 
of the parameter b. The IR fixed points are farest from the singularity if 
b = 2. In this sense this gives the most stable evolution and an optimized 



value for b. Let us note, that the optimized value is usually 6 = 2 [15|, [32 



The evolution of u starts with an irrelevant scaling due to the SG nature of 



the model in the UV [21] . In Fig.[5]one can follow that the evolution turns to 
relevant at about the order of the mass scale y2J. After the relevant scaling 
regime in the crossover the coupling flows to a fixed point marginally. There 
is a certain scale k c , where the marginal scaling appears. This scale can be 
identified with the appearance of the global condensate in the model. The 
value of the IR fixed point depends on z and with line of fixed points in he 
symmetric phase they constitute a continuous line starting from z = I/871", 
u — 1 to z — 1/167T, u = as showed with a thick curve in Fig. [3J The weak 
point of the LPA investigation can be found in the region z > l/8ir. Fig.[6] 
shows some typical evolutions there. Unfortunately the evolution always 
stops since the coupling reaches the singularity condition in Eq. ([8]) . If k 2 < J 
(dotted lines in Fig. [6]) then the coupling scales as u ~ k 2 below the scale 
\plJ therefore the evolution runs into singularity. If k 2 > J the the evolution 
first scales marginally giving universal crossover value (see the solid lines in 
Fig.[6]), but as it reaches the momentum scale V2J the evolution turns to 
relevant giving the flow into singularity again. We note that these scaling laws 
are independent on the scheme, implying that the deep IR regime cannot be 
mapped when z > l/8ir. The results may suggest the splitting of the broken 
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symmetric regime to two different phases as in the massive SG modei 13l. Il5 



and as it was obtained from UV scaling considerations for the LSG model 



331 ] . However this question cannot be answered in the framework of the LPA, 
since the evolutions always run into singularity in there. 

Does the RG method fail to map out the phase structure of the LSG 
model? According to Fig. [3] the answer is affirmative since we have no infor- 
mation on the deep IR scaling when z > I/871". According to fundamental 
considerations the RG flows cannot run into singularity if the effective action 
is not truncated. However we neglected the contributions of the higher har- 
monics, and the gradient expansion in this treatment. The upper harmonics 
can account for the physics of the excitations of higher vorticity, but these 
excitations are usually negligible and we do not expect them to avoid the sin- 



gularity, as in the case of the SG model [2J, [20jj. Furthermore the inclusion 
of the upper harmonics gives only qualitative changes in the flow. Therefore 
one should take into account the next term in the gradient expansion, i.e. 
the wavef unction renormalization. 



3.2. KT transition 

Going beyond the LPA and letting z evolve one obtains the typical phase 
structure of the bi-layer SG model plotted in Fig. [TJ for a certain value of 
Ja = 0.1. The calculations for wavefunction renormalization were performed 
with the choice 6 = 1. Every other parameter value would gave qualitatively 
similar figures. As in the LPA, there is a symmetric phase above l/z\ = 
16tc, and a broken symmetric one below this critical value. The symmetric 
phase contains a line of IR fixed points similarly to the LPA results. It 
implies that the evolution of z is practically negligible in the symmetric 
phase. Furthermore the sensitivity of the effective potential to the UV initial 
values of u also holds. 

The broken symmetric phase drastically changes in comparison with LPA. 
From Fig.[7Jit is clear that there is an attractive IR fixed point of the model 
situated at u — 1, 1/z — and the trajectories tend there. It makes no sense 
to split the broken symmetric phase according to the phase space in Fig. [7] as 
was suggested by the LPA results. Furthermore the figure shows that z — > 00 
in the deep IR regime, which enables to the flow avoid the singularity. In 
Fig. [8] we plotted the flows of the coupling u starting from different UV initial 
values and corresponding to both relations k 2 c > J and k 2 > J. It seems that 
the evolutions run into singularity, however they does not. The inset shows 
that the evolution of z is practically constant and in the vicinity of the 
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singularity they blow up. On one hand it brings the IR value of 1/z — > 0, 
and on the other hand it stops the evolution of the coupling at u = 1. This 
gives the single IR fixed point in the broken symmetric phase, showing that 
opposite to the suggestion of the LPA results, the broken symmetric phase 
is unique. This gives a very serious limitation of the usage of LPA in this 
model. The IR fixed point satisfies the degeneracy condition in Eq. OH]). The 
appearing effective potential in superuniversal, i.e. independent of any UV 
initial parameters. This is a typical property of the spontaneously broken 
symmetric phase. 

The exponent r] deduced from the evolution of z is approximately zero 
in the crossover regime similarly to the SG model. However approaching 
the IR fixed point the evolution of z blows up at a certain scale k c giving 
infinitely large rj. Again, the scale k c is identified with the reciprocal of the 
correlation length £. Its dependence on the reduced temperature t can then 
be easily read off. In case of infinite order phase transition the correlation 
length scales as 

iog£~r". (21) 

In the inset of Fig.[7]we plotted the t dependence of log£, and its slope gives 
v w 0.57 for the exponent. This proves the infinite nature of the phase 
transition of the LSG model. Let us emphasize, that although a KT-type 
phase transition is expected in the LSG model, it has not been proven before. 

4. Summary 

We developed a new strategy in order to identify the type of the phase 
transitions in field theoretical models. It is based on the fact that there is an 
IR fixed point in the broken symmetric phase and the appearing degeneracy 
during the RG flow defines a dynamical momentum scale which reciprocal 
can be identified with the correlation length £, therefore the critical exponent 
v of £ can be determined in an extremely easy way. The value of v may signal 
the corresponding universality class of the model. Our new strategy enables 
us to calculate v in such situations where the RG flow equations has no fixed 
point, too. 

In this article we showed two non-trivial models to demonstrate the power 
of our new method. First we investigated the MSG model, where it is known 
that the inclusion of the mass term in the action changes the phase structure 
considerably as compared to the original, massless one. Our new method 
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showed that the model possesses a second order but not an Ising-type phase 
transition, since we obtained v = 1/2 for the exponent of £. As a second 
example we treated the LSG model. The degeneracy induced scaling strategy 
enabled us to show that the model really has a KT-type phase transition as 
was conjectured earlier. It was also shown that its broken symmetric phase 
is unique. 

The results of o ur p roposed new method may confirm the quantum censor- 
ship conjecture 0, |34j. A proper choice of the suppression and the inclusion 
of the wavefunction renormalization avoids the singularity of the RG flows. 
The scaling of the couplings as the function of k stops in the IR and gives 
a characteristic length to the appearing degenerate vacuum, or condensate, 
while the original momentum scale k reaches zero. 

The developed strategy is capable of describing first order phase transi- 
tions, too. In this case the correlation length does not diverge in the broken 
symmetric phase if we approach the transition point. Thus the type of the 
phase transition can be identified in an easy way according to the scaling 
around the IR fixed point as it was demonstrated in the case of quantum 



Einstein gravity |35J. We note that in this model the existence of the IR 



fixed point has been widely investigated recently [361 ] 
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Figure 1: The scaling of the wavefunction renormalization z according to the two scales k 
(up) and k (down) for the choice b — 2. 
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Figure 2: The reduced temperature dependence of the correlation length for different sup- 
pression schemes, WF refers to the numerical solution including the wavefunction renor- 
malization. 
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Figure 3: The phase structure of the 2-layer SG model is shown in LPA, b = 2. The 
vertical line shows the separatrix at z = 1/167T. The scaling of the coupling u is relevant 
(irrelevant) in the broken symmetric (symmetric) phase, respectively. When z > l/8n we 
have no information about the deep IR scaling. 
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Figure 4: The scheme dependence of the IR fixed point in the case of 1/87T > z > 1/167T 
is presented. The choice 6=1 gives one of the worst, but fastest evolution. The different 
schemes give qualitatively same phase structure. 





Figure 5: The evolution of the coupling u, with the parameter choice z — l/Vl-n and 6 = 2. 
The flows go to the same IR fixed point independently on the initial value of u\ — and J. 
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Figure 6: The evolution of the coupling u, with the parameter choice z = l/An and 6 = 2. 
The flows run into singularity, independently on the on the initial value of Sa = and J. 
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Figure 7: The phase structure of the 2-layer SG model is shown. The trajectories intersect 
since the dimensionless parameter J also evolves. In the inset the scaling of the correlation 
length gives v w 0.57, the straight line with the slope —1/2 is drawn to guide the eye. 
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Figure 8: The evolution of the coupling u with different initial values of uA satisfying 
k 2 c > J (k 2 c > J) and plotted with solid (dashed) line, respectively. In the inset the 
evolution z is drawn. At a certain value of k the flows blow up. 
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